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Abstract 
In  order  to  stuay  the  effects  of  tied  magnetic  lines  in 
ideal  magnetohyarodynamics ,  the  simple  problem  is  treated  of 
the  reflection  of  an  initial  point  disturbance,  i.e.  delta 
function,  froni  a  conducting  wall  in  two-dimensional  magneto- 
hydrodynamics  linearized  about  a  constant  state.   The  lowest 
oraer  magnetic  field  is  taken  normal  to  the  conducting  wall. 
The  conducting  wall  is  the  line  x  =  0  and  the  initial  dis— 
turoance  isatx=x   >0,  y=0.   Tne  solution  separates 
into  three  parts,  the  ordinary  wave  generated  by  a  point 

disturoance  atx=x,y=0,a  reflected  wave  centered  at 

o 

tne  image  point  x  =  -x  ,  y  =  0,  and  a  third  uncentered  wave. 
The  first  reflectea  wave  contains  a  disturbance  similar  to 
the  ordinary  wave  and  a  centered  simple  wave.   The  uncentered 
wave  ana  its  non-self-similar  wave  fronts  are  described.   The 
presence  of  lacunae,  line  singularities,  and  the  nature  of 
the  singularities  at  the  wave  fronts  are  discussed. 


Ill 


I.   Introduction 

It  is  well— known  that  the  boundary  conditions  associated  with 
problems  in  ideal  magnetohydrodynamics  are  very  different  if  the 
boundary  is  a  flux  surface,  n'B  =  0,  or  if  the  magnetic  field 
lines  enter  the  boundary,  n'B  5^  0,  where  n  is  the  normal  to  the 
bounding  surface  and  B  is  the  magnetic  field  vector^   . 
In  the  former  case,  the  boundary  condition  is  n*u  =  0,  where  u  is 
the  flow  velocity  vector,  while  in  the  latter,  usually  referred  to 
as  tied  magnetic  lines,  the  condition  is  u  =  0.   The  case  of  tied 
magnetic  field  lines  is  of  particular  interest  in  stability  analy- 
ses of  various  configurations,  as  it  is  generally  agreed  that 
tied  lines  should  be  much  more  stable  than  untied  lines.   Unfor- 
tunately, there  are  almost  no  concrete  comparisons  of  tied  and 
untied  lines  available.   The  potential  interest  in  tied  magnetic 
lines  leads  us  to  consider  a  much  simpler  problem  in  two  dimen- 
sional magnetohydrodynamic  wave  propagation  of  interest  in  its 
own  right  for  which  explicit  solutions  and  comparisons  are  possi- 
ble.  We  shall  see  that  the  solutions  for  tied  and  untied  lines 
are  totally  different  in  structure. 

We  consider  a  plasma  ±1  the  half— plane,  x  >_  0,  -«  £  y  £  « 
and  we  assume  the  line  x  =  0  to  be  a  perfect  conductor.   We  lin- 
earize the  equations  of  ideal  magnetohydrodynamics  about  a  state 
of  constant  magnetic  field,  plasma  mass  density  and  zero  flow 
velocity.   If  we  take  the  unperturbed  magnetic  field  in  x  direc- 
tion, then  the  field  lines  are  tied  and  the  perturbed  flow  veloc- 
ity vector  vanishes  at  x  =  0.   If  the  unperturbed  magnetic  field 


is  in  the  y  direction,  then  the  field  lines  are  not  tied  and  only 
the  X  component  of  the  perturbed  velocity  vector  vanishes  at 
X  =  0.   We  solve  for  the  Green's  Functions  for  these  problems, 
corresponding  to  giving  delta  function  initial  values  for  the 
perturbations  and  then  following  the  solution  in  time.   The  solu- 
tion for  the  general  initial  value  problem  Is  a  convolution 
Integral  of  the  Green's  Functions  with  the  given  initial  data. 

Some  years  ago  the  Green's  Function  was  given  for  the  case 
of  a  whole  plane  of  plasma  -oo  <_  x  <_  °°,  -co  <^  y  <_  <x>^  and  this  solu- 
tion is  the  basis  of  discussion  for  the  problem  under  consider- 
ation here^  \      Indeed,  the  Green's  Function  for 

the  case  of  untied  field  lines  for  the  half^lane  can  be  trivially 
constructed  by  image  methods.   If  we  place  a  wave  source  at 
x=x   >0,y=Oatt=0  and  a  positive  or  negative  image  source 

atx=-x   <0,y=0,  t=0  and  sum  the  Green's  Functions  from 
o 

these  two  sources,  then  we  can  easily  match  the  boundary  condition 
of  zero  X  component  of  perturbed  velocity.   Thus,  for  untied  mag- 
netic lines  the  Green's  Function  may  be  considered  as  the  sum  of 
the  free  space  Green's  Function  plus  or  minus  Its  reflection  in 
the  plane  x  =  0.   In  particular,  the  only  wave  fronts  present  are 
the  usual  fast  and  slow  wave  fronts  and  their  reflections  in  the 
plane. 

The  Green's  Function  for  the  tied  line  case,  which  we  treat 
in  detail  for  the  rest  of  this  paper,  cannot  be  formed  out  of  the 
free  space  solution  and  images.   We  show,  however,  that  we  can 
split  the  solution  Into- three  waves.   The  first  wave  is  just  the 


free  space  Green's  Function.   The  second  wave  is  centered  on  the 
image  point  to  the  original  disturbance  at  x  =  x  ,  y  =  0;  i.e.  the 
second  wave  is  a  function  of  x  +  x  ,  y  and  t.   Nonetheless,  the 
second  wave  may  not  be  represented  as  a  free  space  wave  with 
source  at  x  =  -x  ,  y  =  0  at  t  =  0.   It  is  composed  of  a  simple 
wave  plus  another  part  whose  wave  fronts  are  just  the  fast  and 
slow  wave  fronts  emanating  from  the  image  point  x  =  x  ,  y  =  0, 
t  =  0.   The  third  wave  is  not  self— similar  in  time  and  has  two 
wave  fronts  and  we  give  figures  showing  these  fronts  developing 
in  time.   We  may  identify  one  of  these  fronts  as  a  wave  which 
leaves  the  wave  source  as  a  fast  wave  hits  the  conducting  wall 
X  =  0  and  reflects  as  a  slow  wave;  the  other  wave  leaves  the  wave 
source  as  a  slow  wave  and  is  reflected  off  the  wall  as  a  fast 
wave . 

The  formal  Fourier  Transform  methods  used  previously  to  de- 
termine the  free  space  Green's  Function  were  not  sufficiently 
sharp  to  Isolate  certain  delta  function  singularities  on  line 
segments,  whose  presence  was  shown  by  exact  plane  wave  methods. 
The  mixed  Fourier- Laplace  Transform  methods  used  here,  which 
automatically  define  distributions  as  boundary  values  of  analytic 
functions,  do  find  these  singularities.   Just  as  for  the  free 
space  case,  the  singularities  ultimately  cancel  out  when  one  spe- 
cializes from  the  general  solution  of  the  fourth  order  partial 
differential  equation  characterizing  these  waves  to  the  special 
type  of   initial  data  peculiar  to  ideal  magnetohydrodynamics . 


One  part  of  the  continuous  spectrum  found  In  magnetohydro— 
dynamic  stability  problems  has  been  associated  with  very  singular 
and  approximately  one— dimensional  wave  propagation  at  certain 
points  of  the  wave  front.   The  behavior  of  the  Green's  Function 
near  this  point  is  analyzed  with  some  care,  and  while  exact  one- 
dimensional  wave  propagation  is  not  found,  the  solution  is  highly 
localized  and  very  singular  there.   Finally,  one  other  interesting 
property  of  the  free  space  Green's  Function  is  the  presence  of 
lacunae,  or  regions  where  the  Green's  Function  is  identically 
zero.   Similar  lacunae  are  found  here  even  though  the  disturbance 
is  reflected  off  a  conducting  wall. 

In  the  next  section  we  formulate  the  precise  mixed  initial 
value— boundary  value  problem  we  intend  to  solve  and  we  represent 
the  solution  in  terms  of  transforms  involving  two  unknown  func- 
tions.  The  third  section  contains  a  detailed  study  of  the  analytic 
properties  of  the  roots  of  dispersion  relation  and  related  func- 
tions.  With  this  analytic  information  we  may  determine  the  two 
unknown  functions  and  complete  the  representation  of  the  transform 
of  the  solution.   All  of  the  above  material  applies  for  general 
solutions  of  the  equations,  but  in  the  remainder  of  the  paper  we 
specialize  to  the  case  of  the  Green's  Function.   In  Section  IV  we 
obtain  implicit  algebraic  representations  of  the  Green's  Function, 
and  the  last  section  explores  various  properties  of  the  solution. 
We  describe  the  development  in  time  of  the  various  wave  fronts, 
and  we  discuss  the  appearance  and  disappearance  of  lacunae  in  the 
various  parts  of  the  solution  and  in  the  whole  solution.   We  fi- 
nally study  the  appearance  of  the  line  singularities  mentioned 
above  and  the  behavior  of  the  solution  near  its  most  singular  points, 


II.   Formulation  of  the  Problem 


We  consider  the  half-^lane  x  >^  0  of  the  x,y  plane  and  we 
treat  the  Lundqulst  equations  of  Ideal  magnetohydrodynamics  lin- 
earized about  a  state  of  constant  density,  pressure  and  magnetic 
field  and  zero  flow  velocity.   We  assume  the  unperturbed  magnetic 
field  to  be  in  the  x  direction.   In  terms  of  the  perturbed  flow 
velocity  vector  (u,v)the  equations  become 


2 
u,.^-a(u   +v   )  =  0 
tt    o   XX    xy 

V. .  -  A^v   -  (a^  +  A^)v   -  a^u   =  0 
tt    o  XX     o    o  yy    o  xy 


(1) 


where  a   is  the  classical  sound  soeed  and  A   is  the  Alfven  speed, 
o  '         o  ^ 

We  note  that  in  terms  of  the  unperturbed  density  p  and  the 
perturbed  density  p 

Po\  +  a^x  =  0         .  (2) 

Since  we  assume  the  boundary  x  =  0  to  be  a  perfect  conductor,  we 
require  that 

u(o,y,t)  =  v(o,y,t)  =  0  (3) 

We  wish  to  solve  the  mixed  initial  value  boundary  value 
problem  for  the  system  (1)  subject  to  the  initial  conditions 

i^(x,y,o)  =  v(x,y,o)  =  0 

u^(x,y,o)  =  f(x,y) 

v^(x,y,o)  =  g(x,y) 

with  f(x,y)  and  g(x,y)  given,  the  boundary  condition  (3),  and  the 
condition  that  no  disturbance  propagates  in  from  infinity.   We 


may  solve  with  more  general  initial  data  by  the  use  of  Duhamel's 
Principle.   In  view  of  (2)  we  see  that  f(x,y)  is  the  x  derivative 
of  some  function.   We  shall  show  below  that  this  property  of 
f(x,y)  will  have  significant  consequences.   After  we  obtain  an 
integral  transform  representation  of  the  solution  of  this  general 
problem,  we  specialize  to  the  case  of  the  fundamental  solution 
characterized  by 

u(x,y,o)  =  v(x,y,o)  =  0 

u^(x,y,o)  =  f6(x  -  x^)<5(y)        f'       ( i|) 

v^(x,y,o)  =  g6(x  -  x^)6(y) 

The  fundamental  solution  is  the  wave  motion  initiated  by  a  point 

disturbance  at  x  =  x  ,  y  =  0  at  time  t  =  0. 

o'  •' 

We  transform  the  system  (1)  as  follows.   We  take  a  Fourier 

transform  in  y  with  transform  variable  k  and  a  Laplace  transform 

in  t  with  transform  variable  p  and  a  Laplace  transform  in  x  with 

transform  variable  s  and  we  define 

2 

^o^x^"'^'^^  ^  a(y,t) 


A^v^(o,y,t)  =  6(y,t) 


so  that 


2   2  2  2 

(p  -a^s  )u(s,k,p)  -  lka^sv(s,k,p)  =  f(s,k)  +  a(k,p)  , 

ika^su(s,k,p)  +  (p^  -  A^s^  +  (a^  +  A^)k^) v(s ,k,p)  =  g(s,k)  +  e(k,p)   , 

where  we  have  used  the  same  symbols  for  functions  and  their  trans- 
forms.  Hence,  we  find 


2  „2  2  -  2,  „2^,  2 


u(s,k,p)  = 


(f(s,k)+a(k,p))  [p'^-A^s  +(a^+A^)k'^)  +  ika^s  (g(s  ,k)  +  6  (k,p) ) 

A(s,k,p) 


(5a) 


ika^s(f(s,k)+a(k,p))  +  (p^-a^s^) (g(s ,k)+B(k,p)) 
v(s,k,p)  =  5 A(s,k,p) '  ^5b) 

where  A(s,k,p)  =  s'^a^A^  -  s^  [a^A^k^+(a^+A^)p^)  +  p"  +  (a^+A^)k^p^ 

(6a) 

or  equivalently  A(s,k,p)  =  p^+(a^+A^)p^ (k^-s^)  -  a^A^Cs^k^-s  )  . 

(6b) 

Now  (5)  and  (6)  do  not  yet  give  the  solution  of  our  problem 

as  (5)  contains  the  unspecified  functions  a(k,p)  and  6(k,p).   In 

order  that  (5)  be  the  transform  of  the  solution  we  require  that 

2   2 
u(x,y,t)  and  v(x,y,t)  must  tend  to  zero  as  x  +y   tends  to  infinity. 

Thus,  u(s,k,p)  and  v(s,k,p)  must  have  no  singularities  in 

Re  s  >  0,  Re  p  >  0.   It  is  clear  from  (5)  that  the  singularities 

of  u(s,k,p)  and  v(s,k,p)  in  Re  s  >  0 ,  Re  p  >  0  come  only  from  the 

zeroes  of  A(s,k,p).   In  the  next  section  we  study  the  zeroes  of 

A(s,k,p)  and  related  functions  and  then  determine  a(k,p)  and  B(k,p) 

so  that  (5)  has  the  appropriate  analytic  structure  to  represent 

the  solution  to  the  original  problem. 


III.   Analytic  Information  and  Determination 
of  the  Transform  of  the  Solution 

We  proceed  to  an  examination  of  the  roots  of  the  dispersion 

relation 

A(s,k,p)  =  0 

considered  as  a  function  of  s  with  k  and  p  considered  parameters. 

We  note  from  (6b)  that  if  s  Is  a  purely  imaginary  root,  then  p  is 

pure  imaginary  as  well.   [This  property  Is  a  direct  consequence 

of  the  hyperboliclty  of  the  system  (1).]   Further,  if  s(k,p)  is  a 

root,  then  so  is  -s(k,p).   In  the  neighborhood  of  p  =  infinity  we 

observe  that  the  four  roots  are  2_  +  o{-] ,   2_  +  o(i),  -  2-  +  0[k] , 

^o  P    o     ^      o     P 

and  -  ?—  +  0 [— ) .   Thus,  if  we  consider  a  root  s(k,p)  as  a  function 

^o     P 
of  the  complex  variable  p,  the  root  is  analytic  in  p  except  at 

branch  points,  has  a  regular  expansion  near  p  infinite,  and  is 

imaginary  only  for  p  imaginary.   Combining  all  this  information 

we  conclude  that  we  may  Introduce  the  four  roots  as 

s-|^(k,p),  S2(k,p),  S2(k,p)  =  -  s^(k,p),  s^(k,p)  =  -  S2(k,p) 

where  i)  the  real  parts  of  p,  s,(k,p),  S2(k,p)  have  the  same  sign, 
11)  s.(k,p),  j  =  1,2,3,4  are  analytic  functions  of  p  except  at  as 

J 

yet  unspecified  branch  points,  and  ill)  near  p  Infinite  s-,(k,p)  and 
s^Ckjp)  have  the  expansions 

^1^^'P)  =  mm  (I    ,A)'   °^^)  ^7a) 

o   o 


'2^^'^^    =max  il    ,A   '   ° (^^    '  ^ ^b ) 


o'  o 


The  case  a  =  A   Is  singular  and  we  must  exclude  It  from  the  dis- 
cussion although  it  is  readily  treated  by  parallel  arguments. 

2 

Since  A(s,k,p)  is  a  quadratic  polynomial  in  s  ,  the  branch 

points  of  s.(k,p)  occur  for  values  of  p  for  which  the  discriminant 
of  the  quadratic  polynomial  vanishes: 

(^'^o^o  -  p'^V^o^'HkW  -  P'(a,-A^)2)  =  0     (8) 
and  the  values  of  p  for  which  s(k,p)  =  0: 

p2(p2  +  k^Ca^+A^))  =  0     .  (9) 

The  branch  points  determined  by  (8) 

ka  A  ka  A 

^        -  a  +A   '     P   -  la  -A  I 
o   o  I  o   o  I 

occur  for  real  values  of  p  so  that  s(k,p)  has  non— vanishing  real 

part  in  the  neighborhoods  of  these  points.   Thus,  if  we  draw  a 

branch  line  between  the  two  branch  points  in  the  half^lane 

Re  p  >  0  and  a  branch  line  between  the  two  branch  points  in  the 

half-plane  Re  p  <  0 ,  then  across  these  lines  s.,  (k,p)  jumps  to 

Sp(k,p)  and  vice  versa  and  s^Ck.p)  jumps  to  S|.(k,p)  and  vice  versa, 

We  turn  next  to  the  branch  points  given  by  (9).   At 

p  =  ^  ik  va  +A  two  of  the  roots  vanish.   For  p  positive  imaginary 

and  IpI  large 

Im  S2(k,p)  <  Im  s-,(k,p)   , 

see  (7).   Thus,  Sp(k,p)  must  vanish  before  s,(k,p)  as  |pl  is  re- 
duced on  the  positive  imaginary  axis.   Hence  we  conclude  that  it 

/   P   p" 
is  SpCkjp)  and  Sh(k,p)  which  vanish  at  p  =  i  ik  /a  +A   .   A  brief 

calculation  shows  that 


(a^+A^)3/\2k)^'^^    .  /n        ro — ?   1/2  . p   1/2 

S2(k,p)  =  ;,  %  ,  .^n   ^     (p-i^^/^ )    ^  o(ip-iic/r;7i) 

(a  +a  A  +A  )^ 
o   o  o   o 

Near  p  =  0  we  obtain  the  expansion  of  the  four  roots  as 

4^+a2 
s  =  1  P  / -o  o°  +  O(p^)  and  s  =  ±  k  +  0(p  ).   Hence,  we  conclude 

o  o 

that  it  can  only  be  s-|^(k,p)  and  s^Ck.p)  which  vanish  at  p  =  0 , 


s^(k,p)  =  pJ-Y~r  ^  °^p^  ^^°^^ 

^0^0 


S2(k,p)  =  k  +  O(p^)  (10b) 

and  that  p  =  0  is  not  a  branch  point.   Thus,  we  may  draw  a  branch 

r~p      2  /~2      2 

line  from  p  =  +  ik  /a  +A   to  p  =  -  ik  /a  +A   and  across  this  line 
'^  vQo  voo 

Sp(k,p)  jumps  to  Sh(k,p)  and  vice  versa.   We  collect  the  above 
information  in  figure  1. 

It  is  useful  to  obtain  a  few  properties  of  the  function 

h(k,p)  =  p^  +  a2s-|^(k,p)s2(k,p)    .  (11) 

Now  h(k,p)  is  well-defined  and  without  branch  points  in  Re  p  >  0 

and  it  may  be  extended  by  analytic  continuation  as  far  as  needed. 

The  points  p  =  ±  ik  /a  +A   are  branch  points  of  h(k,p).   For  lp| 

small  /— P p 

,   /a  +A       p 

h(k,p)  =  pa>  /-5-^  +  0(p'')   ,  (12) 

°  V  a'^A'^ 
o  o 

i.e.  h(k,p)  has  a  first  order  zero  at  p  =  0 .   We  wish  to  determine 
if  h(k,p)  has  any  other  zeroes  in  Re  p  >_  0.   If  h(k,p)  vanishes, 
then  4     4  2      2 


p   -a^s  (k,p)s2(k,p)  =  0   , 


or 


10 


We  already  know  that  h(k,p)  vanishes  at  p  =  0 ,  so  we  must  finally 

check  the  points  ->  o   o   o 

a^k'^Ca^+A^) 


2     o    o   o' 
P   "  ~2~2 


A  -a 
o  o 

?     2 

If  A   >  a  then  the  possible  zero  of  h(k,p)  occurs  for  p  real 
o    o 

and  positive.   In  this  case  the  structure  of  four  roots  is  either 

(s^(k,p),  s^(k,p),  -s^(k,p),  -s^(k,p))  if  s-j^(k,p)  is  complex  or, 

(r-|^(k,p),  r2(k,p),  -r-,^(k,p),  -r2(k,p))  if  r^^Ck.p)  and  r^Ck.p)  are 

real  and  positive.   In  either  case  s-|^(k,p)s2(k,p)  is  real  and 

positive  so  that  h(k,p)  does  not  vanish. 

2     2 
If  A  <  a   ,  then  the  possible  zero  of  h(k,p)  occurs  for  p  pure 

2 
imaginary.   Since  p   is  there  real  and  negative  and  s^(k,p)  is  pure 

imaginary  h(k,p)  can  vanish  only  if  Sp(k,p)  is  pure  imaginary  and 

Im  s,(k,p)  and  Im  s„(k,p)  are  of  the  opposite  sign.   Wherever 

Sp(k,p)  is  pure  imaginary  Im  s,(k,p)  and  Im  S2(k,p)  have  the  same  sign, 

Hence,  h(k,p)  vanishes  in  Re  p  >^  0  only  at  p  =  0. 

We  are  now  prepared  to  return  to  the  specification  of  a(k,p) 
and  g(k,p)  in  (5)  so  as  to  complete  the  solution  of  our  problem. 
We  recall  that  u(s,k,p)  and  v(s,k,p)  should  be  analytic  in  Re  s  >  0 
and  Re  p  >  0,  but  that  their  representations  in  (5)  and  (6)  indi- 
cate that  they  are  singular  at  the  zeroes  of  A(s,k,p).   We  know 
that  A(s,k,p),  considered  as  a  function  of  s,  has  exactly  two  zeroes 
in  Re  s  >  0,  s-]^(k,p)  and  S2(k,p).   We  select  a(k,p)  and  3(k,p)  so 
that  u(s,k,p)  and  v(s,k,p)  are  regular  at  s,(k,p)  and  S2(k,p), 
i.e.  we  require 
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(f(Sj.(k,p),k)+a(k,p))(p^-A^s2(k,p)  +  (a2+A2)k2) 

+  ika^s.(k,p)(g(s.(k,p),k)  +  3(k,p))  =  0  (13) 

O  J  J 

J  =  1,2 
and 
Ika^s  (k,p)(f(s.(k,p),k)+a(k,p)) 

+  (p2-a2s?(k,p))(g(s.(k,p),k)+6(k,p))  =  0   ,     J  =  1,2  .(14) 

We  see  that  we  have  four  equations  and  two  unknown  functions. 
However,  a  simple  computation  shows  that  since  A  (s-|^(k,p)  ,k,p)  =  0 
the  two  equations  (13)  j  =  1  and  (l^J)  j  =  1  are  linearly  depen- 
dent.  Similarly  (13)  j  =  2  and  (l4)  j  =  2  are  also  linearly  de- 
pendent.  Thus,  we  have  eliminated  the  undesired  singularities  if 
we  satisfy  (13)  or  (l4).   Choosing  (1^) ,   we  find 

-(^'P)  =  is,(k,p)-sjk,p)jh(k,p)  (3i(^,P)f(3,(k,p),k)[p^-a^s^(k,p)) 
-  S2(k,p)f(s2(k,p),k) (p^-a^sj(k,p)) 


2  2 

+  2-^  fg(s.(k,p),k)  -  g(s,(k,p),k))  ^  (15a) 

o 


(g(s2(k,p),k)  -  g(s^(k,p),k))  I 


e^^'P)  =  is2(k,p)-s,(k,p)Jh(k,p)  {lksi(k,p)s2(k,p) 
.  [f(s^(k,p),k)-f(s2(k,p),k)) 
+  s^(k,p)(p^-a^S2(k,p))g(s2(k,p),k} 
-  S2(k,p)[p^-a^s^(k,p))g(s^(k,p),k)  V  (15b) 
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We  may  readily  verify  that  a(k,p)  and  g(k,p)  are  analytic  in 
Re  p  >  0  and  have  no  singularities  at  the  branch  points  of  s.(k,p) 

J 

Further,  they  may  be  extended  onto  Re  p  =  0  with  no  singularity 
at  p  =  0.   Thus,  with  a(k,p)  and  6(k,p)  as  given  by  (15),  we  con- 
clude that  (5)  represents  the  transform  of  the  solution  to  our 
problem. 
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IV.   A  Representation  of  the  Fundamental  Solution 

We  now  turn  to  the  fundamental  solution  characterized  by  ( ^l ) 

so  that  _3jj 

f(s,k)  =  fe   ° 

-sx 
g(s,k)  =  ge 

Since  the  representations  are  rather  cumbersome,  we  carry  through 

the  details  only  for  u(x,y,t).   Exactly  the  same  procedures  apply 

for  v(x,y,t)  as  well,  although  v(x,y,t)  is  simpler  in  one  respect 

which  we  note  at  the  appropriate  point.   If  we  Invert  the  s 

transform,  we  find  different  forms  for  x  <  x^  and  x  >  x^ .   For 

X  >  x  we  conclude 
o 

H  sj(k,p)x  -s.(k,p)x 

-^-'^'P)  =  .^3  A^iL(k,p),k,pJ  [^'^     '  °  ^  "^^'P)) 

.  (p2-A2s2(k,p)  +  k^Ca^+A^)} 

^  -s.(k,p)x  1 

+  ika^s_.(k,p)[ge   J      °  +  B(k,p))  >    ,       (l6a) 

while  for  x  <  x   we  have 

-(-'^'p)  =  J^  A3i:.(k,p),k,pj  [aiKP)[p'-^y^i^,p)  -  ^'^-K^) 

+   ika^sj(k,p)e(k,p)  [  .  (I6b) 

We  could  verify  explicitly  from  (16)  that  u(x,k,p)  tends  to  zero 
as  IpI  -»-<»,  Re  p  >  0,  as  a  Laplace  transform  must. 

We  may  combine  (l6a)  and  (l6b)  into  one  formula  by  the  use  of 
(13)  and  the  relations  s-|^(k,p)+ST  (k,p)  =  0  =  S2  (k,p  )+S2^  (k,p) 
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^        Sj(k,p) |x-x^|   ^ 
J  -  J   s   J  I, 


+  ika^s.(k,p)g  sgn(x-x^) 
s  .  ( k ,  p  )  X 


^3  A3[s.(k,p),k,pJ  {  «(K>P)(p'-A^j'(k,p).k^(a^A^)) 


+  ikaS  (k,p)B(k,p) 


(17) 


-s  (k,p) 
Since  a(k,p)  and  6(k,p)  are  linear  combinations  of  e 
-S2(k,p)x^  s-,(k,p)x^      S2^(k,p)x^ 

and  e  or  equivalently  e  and  e  we  see 

that  u(x,k,p)  is  a  linear  combination  of  a  function  of  x-x  ,  a 

function  of  x+x  ,  and  a  more  general  function  of  x  and  x  .   In 

(17)  the  function  of  x-x   is  given  by   the  first  summation  and  the 

other  two  functions  are  given  by  the  second  summation.   We  set 

u(x,y,t)  =  u-|^(x-x^,y  ,t)  +  U2(x+x^,y,t)  +  U2(x,x^,y,t)   , 

where  the  transform  of  u-,  is  given  by  the  first  summation  in  (17) 
and  the  sum  of  the  transforms  of  u^  and  Uo  is  given  by  the  second 
summation  in  (17)-   We  treat  the  terms  separately  and  start  with 
the  simpler  u, (x-x  ,y,t). 

The  representation  of  u, (x-x  ,y,t)  is 


u-|^(x-x^,y,t)  =  Re 


277  1 


«>     0  +  1°° 

dk 
0   a-l^o 


^    Sj(k,p)|x-x^ 


,    iky  pt   V   e 

'P  '        '        ^3  A3is.(k,p),k,pJ 


(18) 


<    f(p^-A^sJ(k,p)+k^(a^+A^))  +  ika^Sj(k,p)g  sgn(x-x^) 
where  a  >  0.   We  now  deform  the  contour  of  integration  in  the 
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K7 

p  plane  to  CI  or  C2  [see  Figure  2]  as  t  -  /  ^  g"  I ^~^o '  ^^  greater 

o  o 
than  or  less  than  zero.   We  readily  verify  that  on  these  contours 

Re  (pt+s . (k,p) |x-x  I)  10,  j  =  3,^-   If  we  denote  the  integrand  in 

(18)    as    I,   we    find 

u-L(x-x^,y,t)    =   Re   — 


2TT^i 


dk 


dp    I 


0      CI   or  C2 


27r 


dk   e 


Iky 


/a^A^  /  /a^A^ 


a  A' 
o   o 


a   A' 
o   o 


In  the  first  integral  we  Introduce  a  new  variable  p   =  p/k  and  we 
find  that  we  may  perform  the  k  integration  by  means  of  the  fol- 
lowing distribution  identity 

ika 


dk  e- 


TT(5(a)  +  i/a 


where  i/a  is  understood  to  be  in  a  principal  value  sense: 


u-,^(x-x^,y,t) 


=  Re 


2.h 


CI  or  C2 


4  r  f [p^-A^s^(l,p)+(a^+A^)+ia^s.(l,p)g  sgnCx-x^)) 


s.(p) 
•  (TT6(y+  I  t  +  -^^—    |x-x^|))  + 


s.(p) 

V  +  ?-  +  -^. Ix-x 

"     1      1    '    o 


o  o  \  o  o 
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The  delta  function  under  the  Integral  sign  does  not  contribute  to 

the  integral,  for  either  it  has  non— zero  argument  or  where  its 

argument  is  zero,  the  contribution  is  pure  imaginary.   We  deform 

the  p  integration  to  a  contour  just  to  the  right  of  the  imaginary 

s  .  (p)        _-, 
axis.   We  recall  that  (y  +  ?■  +  -^ —  '  ^"■'^n  ' -^    "^^  defined  in  a 

principal  value  sense,  so  we  must  add  in  one— half  the  residue  at 
the  points  at  which  it  has  a  pole.   As  with  the  delta  function, 
the  residues  yield  a  purely  imaginary  term,  which  does  not  con- 
tribute to  u, .   Thus,  we  obtain 


u-,^(x-x^,y,t) 


o+i-    4   f(p2-A2s2(i,p)+(a^+A2))  +  ia^gs . (1 ,p)sgn(x-x^ ) 


Re 


2tt 


0-1° 


dp  I 
.        0  =  3 


o  J 


o"  J 
s,(l,p) 


A^(s.(l,p),l,p)(y  +  £  +  -^ l^-^oD 


o  o 


t  -   x-x 


a  +A 
o o 

a  A 
o  o 


where  a  >  0  and  there  are  no  singularities  in  0  <  Re  p  <_  a.   Fi- 
nally we  note  that 

s.(l,p) 
Im  (y  +  f  +  -L |x-x^|)  <  0 

for  Re  p  >  0  and  j  =  1,2.   Hence  we  nay  add  into  the  summation 
the  terms  with  j  =  1,2,  after  observing  that  the  branch  line  in 
Re  p  >  0  disappears  when  we  combine  the  two  terms  with  j  =  1  and 
j  =  2.   We  find 


u-|^(x-x^,y,t)  =  Re  — ^ 

2Ti 


°^^°°  ,  f(P^-AoC^(p)+a^A^)  +  ia^gap)sgn(x-x^) 
'^^  A[?(p),l,pJ|x-x^| 


0_loo 


/  2   2 
o      o 


a'^+A'^ 
o   o 


..y_j_..<„„  ,_y^,_^, 


o  o 


(20) 
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where 

The  representation  (20)  is  essentially  the  same  as  found  earlier, 
although  expressed  in  other  variables  [Weitzner,  196I].   Thus, 
(20)  is  just  the  fundamental  solution  for  a  point  disturbance  at 
X  =  X  ,  y  =  0,  t  =  0.   As  was  shown  earlier  (20)  admits  an  alge- 
braic representation  in  terms  of  the  zeroes  in  p  of  A(c(p),l,p), 
a  fourth  degree  polynomial  in  p. 

The  presence  of  the  delta  function  in  (20)  is  surprising  on 
two  counts.   The  transform  methods  used  earlier  were  not  sharp 
enough  to  pick  out  the  delta  functions,  whose  presence  was  shown 
by  plane  wave  methods.   The  techniques  used  here  have  been  able 
to  isolate  the  delta  functions.   The  other  surprise  is  the  presence 
of  any  delta  function.   As  was  shown  earlier,  there  should  be  no 
delta  functions  for  the  special  case  of  magnetohydrodynamic  wave 
motion,  as  opposed  to  more  general  wave  motion  satisfying  the  same 
fourth  order  partial  differential  equation.   With  regard  to  the 
first  point  we  recall  that  one  way  to  define  a  distribution  is  as 
the  limiting  value  of  an  analytic  function.   By  the  use  of  Laplace 
transforms,  with  natural  analytic  properties,  we  might  readily 
extend  the  fundamental  solution  into  the  complex  x,y  plane  and 
thereby  define  the  distributions  properly.   The  Fourier  transforms 
■used  previously,  confuse  analytic  properties  in  x  and  y  and  do 
not  admit  of  direct  extensions. 

With  regard  to  the  second  point,  we  can  readily  show  that 
the  singularity  on  the  full  line  segment  of  support  of  the  delta 


p7? 

function,  y  =  0,  t  >  J   p   2°    I ^"^n ' '  ^^  "°^  present.   We  recall 

o  o 
from  (2)  that  for  any  magnetohydrodynamics  problem  f(x,y)  =    ^^'^ 

Thus,  if  we  take  the  convolution  of  the  delta  function  part  of 

the  solution  with  f(x,y)  we  find 


dx  dy  f(x-x  ,y-y  )H|t  -  j    2    2°    '^'"^oM  ^^^    ^ 

^0^0 


\  O    O      /       \  0    0/ 

Thus,  for  a  magnetohydrodynamics  problem  only  the  end  points  of 
the  above  line  segment  are  in  the  support  of  the  distribution. 

We  analyze  the  behavior  of  the  fundamental  solution  near  the 

r"2   2 
a  +A 


points  y  =  0,  t  -  /  ?,      °    |  x-x  |  =  0  in  greater  detail  in  the  next 

^a-A'^      ° 
o  o 

section  and  we  show  that  if  one  solves  a  magnetohydrodynamics 
problem  for  which  the  initial  data  f(x,y)  =   '^^^'^  ,  then  the 
delta  function  exhibited  here  explicitly  is  canceled  by  singular- 
ities present  in  the  remainder  of  the  fundamental  solution.   Sim- 
ilarly, no  delta  functions  are  present  in  the  part  of  the  funda- 
mental solution  associated  with  the  data  g(x,y).   These  complex- 
ities are  all  absent  in  the  function  v(x,y,t)  for  which  the  anal- 
ysis is  not  carried  out. 

We  return  to  (17)  and  follow  the  same  procedure  to  study 
Up(x+x  ,y,t)  and  Uo(x,x  ,y,t).   We  introduce  the  notation 
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-s, (k,p)x 


a(k,p)  =  0LAk,p)e      -^  ""  +  a^(k,p)e 

-s-^(k,p)x^ 


-s^{k,p)x^ 


B(k,p)  =  62(k,p)e 


+  B^(k,p)e 


-S2(k,p)x^ 


where 


a^Ckjp) 


f(p^-a^sg(k,p))s3(k,p)  -  ikp^a^g/A^ 
[S3(k,p)-s^(k,p)jh(k,p) 

f(p2-a^s2(k,p))s^(k,p)  -  Ikp^a^g/A^ 
"il^^'P^  =  ls,(k,p)-s,(k,p)Jh(k,p) 


and 


B3(k,p)  = 


e4(k,p) 

We  now  find 


o  0       0    0 

ika^S3(k,p)sn(k,p)f  +  (p  -a^s^Ck.p))  s^^  (k,p)g 

(S3(k,p)-s^(k,p)jh(k,p) 

O  0         0     0 

ika^S3(k,p)s^(k,p)f  +  (p  -a^s^ (k,p)) S3(k,p)g 


[s^(k,p)-S3(k,p))h(k,p) 


(22a) 


,  (22b) 


(23a) 


(23b) 


(24a) 


(24b) 


U5(x+x  ,y,t)  +  u,(x,x  ,y,t)  =  Re  — p- 


iky  , 
e  •^dy 


a+1'^ 


a-i" 


ePtdp 


\X   I  rs.ik,p),k!p]  [aj(k,p)(p2-A2s2(lc,p).k2(aX)>ll<ah(^.P)Pj(^'P)^ 

s-(k,p)x+S2|(k,p)x^ 
^  "    A  (S3(k,p),k,p)    [«4(k.P)(p'-A^s2(k,p).k2(^2,;,2)^,,j^^23^(j^^p)3^(,^p)3 


2  „2  2, 


2,  2.  „2. 


s^(k,p)x+S3(k,p)x^ 
'  "   Ajs^(k,p),k,p]    [«3(k,P)(p'-A^s^(k,p).k^(aX))^ll^<-4('^'P)e3^^'P^^^ 

(25) 
We  see  that  Up(x+x  ,y,t)  is  given  by  the  summation  in  (25)  and  that 
u-,(x,x  ,y,t)  is  given  by  the  last  two  terms.   We  note  that  the 
transforms  of  Up  and  u-,  have  poles  at  the  branch  points  of  s(k,p) 
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in  Re  p  >  0  although  the  sum  of  the  transforms  is  regular  at  these 
points.   Thus,  we  do  not  separate  the  two  sets  of  terms. 

If  we  denote  the  integrand  of  (25)  by  I,  and  we  move  the  path 

of  integration  in  the  p  plane  to  CI  or  C2  as  t  -  /  S  p°  (^"•"^o^  ^  '^ 
'  2,„2  ''  alA, 


/a"+A 
or  t  -  /  °  °    (x+x  )  <  0,  we  find 

o  o 


o  o 


u,(x+x  ,y,t)u  (x,x  ,y,t)  =  Re  — x- 


dk 
0   CI  or  C2 


dp  I 


o  o    \     o   o 


We  may  now  perform  the  k  integration  after  introducing  p  =  p  k. 
Vie   recall  that  singular  denomenators  are  to  be  taken  in  a  princi- 
pal value  sense.   We  next  move  the  path  of  integration  just  to  the 
right  of  the  real  axis  and  we  find  the  representation 
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a+i°° 


u„(x+x^,y,t)+u_(x,x  ,y,t)  =  Re 


dp 


2tt' 


a-i° 


^      a^(l.p)(p^-A^s_^(l,p)+a^+A^)  +  B_.  ( 1  ,P  )la^s  .  (1  ,p) 

,1^  s.(l,p) 

J  =  3    Ajsj(l,p),l,p)(y  +  ^  t  +  -L_ (x+x^)) 

a^(l,p)(p^-A^s^(l,p)+(a^+Ag))  +  la^Sg(l,p)B^(l,p) 

■^                    ~   s  (l,p)x+s.(l,p)x 
A^[s  (l,p),l,p](y  +  I  t  +  -3 ^ ^  ) 


2  „2_2 


2.  „2> 


a,(l,p)(p'^-A^st(l,p)+a^+A^)  +  ia  s  .  (1  ,p)  B  (1  ,p) 


o   o' 


'*'  IT    S2.(l,p)x+s,(l,p)x 

A3(s^(l,p),l,p)(y  .21  ,  Ji _J 03 


i  7a^+A 


+  Re  ^r 

TTl 


O   O 
dp 


2  „2  2 


2.  „2> 


.  2 


a^ ( 1 ,P )  (P  -AqS  2 ( 1 , P ) +aQ+A^) +B2(l,p)iaQS^(l,p 


-1  /a^+A^ 
o   o 


Ag(S2(l,p)  ,l,p] 


s,(p)  as,(l,p) 


1-   s,vp ;  c 


a  A' 
o  o 


a  A' 
o  o 


We  have  freely  used  the  reality  properties  of  the  various  functions 
to  simplify  the  expressions  and  obtain  (26).  The  first  integral  is 
taken  on  the  path  Re  p  =  o  >  0,  where  there  are  no  singularities  of 
the  integrand  in  a  >^  Re  p  >  0. 

In  (26),  u^(x,x  ,y,t)  is  given  by  the  second  and  third  terms 
of  the  first  integral;  the  remainder  of  (26)  gives  U2 (x+x^  ,y ,t ) . 
As  (26)  is  written,  there  is  no  singularity  of  the  integrand  on 
the  branch  line  of  s.(l,p)  in  Re  p  >  0  although  there  are  poles 

J 
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at  the  branch  points  at  the  ends  of  the  branch  lines.   The  poles 

in  the  integrands  of  Up  and  u^  cancel  each  other.   We  may  then 

consider  Up  and  u-,  separately,  provided  we  ignore  the  contribution 

to  the  contour  integral  from  the  spurious  poles.   From  (26), 

UpCx+x  ,y,t)  is  the  sum  of  three  terms,  a  contour  integral,  the 

Integral  of  a  delta  function,  and  a  delta  function  alone.   The 

lone  delta  function  is  similar  to  that  found  in  u-,  (  |  x-x  |,y,t) 

and  in  the  next  section  we  show  that  for  magnetohydrodynamics 

problems  it  is  canceled  by  other  singularities  in  Up .   We  may 

readily  evaluate  the  integral  of  the  delta  function,  the  second 

term,  and  we  find  that  the  structure  of  this  term  is 

ds-,(l,p) 
(b(p)H(t  +    (x+x^)    ^p      ) 

ds^(l,p) 

ft  +  (x+x  )  — i ] 

*■        o    dp      -" 

x+x 
where  p  is  a  function  of  — —  and  ^  determined  by 

x+x         x+x  x+x 


and  the  conditions  that  p  be  pure  imaginary  and 


-i  Va  +A   <  Im  <  i  /a  +A 
^  o   o  —    —   vy  o 


ds  (l,p)  _,              ds-,(l,p) 
Up  to  the  factor  (t  +  (x+x^)  — ^ }    H[t  +  (x+x^)  — ^ ], 

this  term  is  constant  on  a  plane 

ly  +  pt  +  {x+x^)s^{l,p)    =  0   , 

where  p  is  fixed.   Thus  we  identify  this  term  as  a  simple  wave. 
The  remaining  part  of  Up(x+x  ,y,t)  is  the  contour  integral.   The 
singularities  of  this  integral  in  Re  p  >  0  occur  at 
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ly  +  pt  +  s  .  (l,p)  (x+Xq)  =  0,  j  =  3,4.   Just  as  with  u-j^,  for  Re  p  >  0 , 
Re  (ly  +  pt  +  s.(l,p)(x+x  ))  >  0,  j  =  1,2,  so  that  the  singular!- 
ties  of  the  integrand  of  the  contour  integral  contribution  to  u^ 
occur  at 

A(ifif  .i.P)  -  0   . 

Hence,  we  may  express  this  part  of  U2  explicitly  and  algebraically 
in  terms  of  residues  of  the  Integrand  at  the  roots  of  this  fourth 
order  equation. 

As  was  just  noted,  Uo(x,x^,y,t)  is  given  by  a  contour  integral 

with  the  second  and  third  terms  of  the  contour  Integral  in  (26). 

The  poles  of  this  integrand  that  contribute  to  U2+U2  occur  at  the 

roots  in  Re  p  >  0  of 

ly  +  pt  +  S2(l,p)x  +  s^(l,p)x^  =  0 

and 

ly  +  pt  +  s^(l,p)x  +  s^(l,p)x^  =  0    , 

and  u-j  is  then  represented  in  terms  of  the  residues  of  the  inte- 
grands at  these  values  of  p.  We  note  that  these  roots  are  among 
the  zeroes  of  the  eighth  degree  polynomial 

^    Ui^I^£^,l,^),    A(^,l,p))/[A(ig£^,l,p)  A[i^,l,p)] 
^     ^o  o  o  o 

where  R  ff  (tt)  ,e(TT)l  is  the  eliminant  with  respect  to  ir  of  the  two 

polynomials  f(Tr)  and  gCir).   Simpler  explicit  representations  of 

the  polynomial  are  possible.   In  the  next  section  we  explore  some 

properties  of  the  fundamental  solution  just  obtained. 
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V.   Some  Properties  of  the  Fundamental  Solution 

In  this  section  we  explore  some  of  the  properties  of  the 
fundamental  solution.   The  surfaces  on  which  the  solution  is  sin- 
gular are  the  wave  fronts  and  we  obtain  the  equations  for  them. 
We  can  determine  various  limiting  points  on  the  fronts  and  we  ob- 
tain the  geometrical  connections  amongst  the  various  surfaces. 
We  also  relate  the  simple  wave  to  these  wave  fronts.   Figures  ob- 
tained by  direct  computation  show  the  development  of  the  fronts  in 
time.   V/e  also  observe  the  presence  of  lacunae  in  the  solution. 
Finally,  we  study  in  some  detail  the  behavior  of  the  fundamental 
solution  in  the  neighborhood  of  the  very  singular  points  y  =  0, 

a  A 
o  o 

If  we  consider  representations  (19)  and  (26)  for  parts  of  the 

fundamental  solution,  we  may  readily  verify  that  the  solution  is 

smooth  in  x,  y,  and  t  except  at  points  where  the  zeroes  in  p  of 

iy  +  pt  +  Sj(1,p)|x-Xq| ,     j  =  3,^  (27) 

or         iy  +  pt  +  s.(l,p) |x+x^| ,     j  =  3,^  (28) 


or 


iy  +  pt  +  S2(l,p)x  +  s^(l,p)x^  (29) 


or         iy  +  pt  +  3^(1, p)x^  +  s^(l,p)x  (30) 

move  from  the  half— plane  Re  p  >  0  onto  the  imaginary  axis. 
Re  p  =  0 .   If  we  denote  one  of  the  above  functions  by  h(x,y,t,p) 
then  the  conditions  that  a  root  in  p  move  onto  the  imaginary  axis 
are  just 
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h(x,y,t,p)  =  0 


||(x,y,t,p)  =  0 


(31) 


where  p  is  pure  imaginary.   At  any  given  time  (31)  is  a  parametric 
representation  of  a  curve,  a  wave  front.   Thus,  if  h(x,y,t,p)  in 
(31)  is  replaced  successively  by  the  functions  in  (27)  through 
(30),  we  obtain  the  several  wave  fronts  in  this  problem.   This 
characterization  of  the  wave  fronts  is  identical  with  the  usual 
geometrical  optics  description  of  wave  fronts  of  hyperbolic  par- 
tial differential  equations    .   We  note 

that  for  (27)  and  (28)  with  j  =  3,  P  ranges  over  the  entire  imag- 
inary axis,  while  for  (27)  and  (28)  with  j    =   ^   and  for  (29)  and 


(30)  p  is  Imaginary  and  |  Im  p  |  >_  /a  +A 


If  we  consider  the  function  u,(x-x  ,y,t),  whose  wave  fronts 
are  determined  by  (27)  we  find  that  they  appear  as  in  figure  3, 
where  we  have,  for  convenience,  set  x   =  0.   The  outer  curve  is 
called  the  fast  locus  and  is  given  by  (27)  with  j  =  4,  while  the 

inner  curves  are  called  the  slow  locus  and  are  given  by  (27)  with 

/  2„2 
a   A 

J  =  3.   We  refer  to  the  points  y  =  0,  x  =  ±t  /-§— %  as  the  cusps, 

•'a  +A'^ 
o   o 

although  other  cusps  are  obviously  present.   We  use  these  terms 
to  describe  the  other  wave  fronts  as  well.   Inside  the  inner  curves 
we  have  shown  previously  that  the  fundamental  solution  is  zero. 
Thus,  these  regions  are  lacunae,  or  gaps,  in  the  solution. 

Now  U2(x+x  ,y,t)  as  given  by  (26)  Is  the  sum  of  three  terms, 
a  contour  integral,  a  simple  wave  coming  from  the  integral  of  the 
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delta  function,  and  the  delta  function.   The  part  of  u„(x+x  ,y,t) 

excluding  the  simple  wave  has  essentially  the  same  wave  fronts  as 

u-,(x-x  ,y,t)  except  that  the  former  fronts  are  centered  on 

X  =  -X  ,  y  =  0  while  the  latter  fronts  are  centered  on  x  =  +x  , 
o  o 

y  =  0.   We  may  characterize  the  fronts  of  Up(x+x  ,y,t)  as  the 

front  of  a  fast  wave  which  leaves  x=x  ,y=Oatt=0,  reflects 

off  the  conducting  wall  as  a  fast  front,  (28)  with  j  =  ^,  and  the 

front  of  a  slow  wave  which  leaves  x=x  ,y=Oatt=0,  reflects 

o 

off  the  conducting  wall  as  a  slow  front,  (28)  with  j  =  3-   It  is 
clear  from  (20)  and  (26)  that  while  this  part  of  the  solution  is 
is  centered  at  the  image  point, x  =  -x  ,  y  =  0,  to  the  original 
disturbance,  x  =  +x  ,  y  =  0 ,  it  is  not  merely  the  reflection  of 
u, (|x-x  |,y,t),  but  a  more  complicated  function. 

We  may  readily  determine  that  in  the  same  light  (29)  repre- 
sents a  wave  front  which  leaves  x=x  ,y=0att=0asa  fast 
wave,  and  is  reflected  off  the  conducting  wall  as  a  slow  wave. 
Similarly  (30)  represents  a  wave  which  leaves  x  =  x  ,  y  =  0  at 
t  =  0  as  a  slow  wave  and  reflects  off  the  conducting  wall  as  a 
fast  wave.   These  are  the  wave  fronts  of  Uo(x,x  ,y,t). 

We  may  readily  determine  the  position  of  certain  limiting 
points.   From  (10a)  we  conclude  that  as  p   tends  to  zero,  (27) 
with  j  =  3  determines  the  points  on  the  slow  wave  front  of  u-,  as 

y  =  0,    |x-x  I  =  /-f-S"  ^> 

^  a  +A 
o   o 

i.e.  the  cusps,  while  (28)  with  J  =  3  yields  the  point  on  the  slow 

locus  of  u^ 
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/a  A 
y  =  0,      x+x   =  /-§-^  t, 

o   o 
again  the  cusp.   At  the  other  limit,  as  p  tends  to  infinity  we 
obtain  from  (7a)  the  points  on  the  slow  wave  front  of  u, 

y  =  0,      I x-x  I  =  t  min  (a  ,A  ) 

J  y  '     O  '  O  '  O 

and  the  points  on  the  slow  wave  front  of  Up 

y  =  0,     x+x   =  t  min  (a  ,A  ) 
''     '        o  o '  o 

Some  of  these  points  may  lie  in  x  <  0  for  t  sufficiently  small 

and  are  irrelevant. 

For  the  remaining  wave  fronts  the  range  of  admissible  values 


y2   2 
a  +A  .   Thus,  it  is  not 


yp   p 
a  +A 

we  see  that  | -; 1  ->  »  so  that  no  point  in  x  >_  0  is  defined 

by  (29)  while  for  (30),  i.e.  for  the  slow  converted  to  fast  front 


of  u^,  we  obtain  the  points 


,  5      2    2      5" 
'a  ^a'^A'^+A 


t  J4<  -  \  r  I2  °]  . 


a   A 
o   o 


provided 


ds    (l,p)  /-^ — p 

+   X      — i >    0      at      p   =   ±i    /a   +A        . 

o        dp  —  ^  >/    o      o 


dp 
If  the  inequality  fails  then  no  point  is  determined  by  (31)  as 


p  ->  ±i  Ja  +A  .   We  may  take  the  limit  as  p  tends  to  infinity  and 
we  find  that  points  on  the  fast  locus  of  u-,  are 

y  =  0,      l^-^o'  ^  ^  "^^^  ^^o'^o^ 
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an 


d  a  point  on  the  fast  locus  of  u„  Is 


y  =  0      x+x   =  t  max  (a  ,A  )    , 
and  points  on  the  locus  of  u-,  are 

y  =  0      X  max  (^,  j^]    +   x^  min  (— ,  ^)  =  t 

O     O  0     0 

and      y  =  0      x  min  [i-,  i-]  +  x  max  (— ,  ^)  =  t 

O    o  o    o 

Clearly  the  former  is  a  fast  wave  converted  to  a  slow  wave  at  the 
conductor  and  the  latter  is  a  slow  wave  converted  to  a  fast  wave 
at  the  conductor. 

We  see  that  certain  wave  fronts  coincide  at  x  =  0 .   The  fast 
fronts  of  U-,  ,  u^  and  the  fast  converted  to  slow  front  of  u^  coin- 
cide at  X  =  0.   Similarly,  those  of  the  following  fronts  that 
touch  X  =  0  coincide,  the  slow  fronts  of  u, ,  Up,  and  the  slow 
converted  to  fast  front  of  u^-   We  might  show  the  possibility  of 
the  appearance  and  disappearance  of  cusps  in  the  wave  fronts  of 
Ut,  but  we  leave  the  analysis  to  the  reader  and  illustrate  the 
phenomenon  by  computed  figures. 

We  next  consider  the  geometry  of  the  simple  wave  characterized 

by 

iy  +  pt  +  s-,(l,p)(x+x  )  =  0 

ds  (l,p) 

and  t  +  — ^- (x+x^)  >  0 

dp         o  — 

2   2" 
a  +A  .   We  note  that  the  simple 

wave  disappears  at  points  at  which  the  inequality  becomes  equality. 

Thus,  the  simple  wave  ends  on  a  wave  front,  specifically  the  slow 


29 


wave  front  of  Up .   Further,  as  x  ->  +"  the  inequality  is  violated, 
so  that  the  simple  wave  lies  to  the  left  of  the  slow  wave  front 
of  Up.   Recalling  that  the  simple  wave  is  essentially  constant  on 
the  lines  given  by  the  equality,  we  see  that  the  line  of  largest 


•  /r2,  „2 

=    X  Ja 


slope  has  p  =  i  Ja   +A   and  the  intercept  of  this  line  with  x  =  0 


is  the  intersection  with  x  =  0  of  the  slow  converted  to  fast 
locus  of  u-,.   Further,  the  slopes  range  in  absolute  value  from 
zero  to  this  maximum. 

We  exhibit  in  the  following  sequences  of  figures  the  behavior 
of  the  wave  fronts.   We  do  not  plot  the  fronts  of  u,  as  they  would 
complicate  the  pictures.   We  illustrate  the  case  in  which  the 
ratio  of  the  sound  speeds  is  .8.   The  figures  are  drawn  to  differ- 
ent scales.   The  position  of  the  initial  singularity  x  =  x  ,  y  =  0 

is  marked  with  an  x  on  each  figure  and  the  time  is  given  in  units 

x 

of     °   .   In  figure  h    the  fast  locus  of  u„  and  the  fast  converted 

/a  A  ^ 

o  o 

to  slow  locus  of  u^  appear.   Their  intercepts  with  x  =  0  coincide 
and  the  outer  locus  is  the  fast  locus  of  Up.   In  figure  5  the  slow 
locus  of  Up  and  slow  converted  to  fast  locus  of  u^  appear,  al- 
though it  is  difficult  to  separate  the  curves.   In  figure  6  the 
slow  locus  of  Up  and  the  slow  converted  to  fast  locus  of  u^  are 
easily  separated,  the  latter  extending  to  larger  values  and  not 
intersecting  x  =  0.   Note  that  the  cusp  of  u,  has  not  yet  touched 
X  =  0.   After  the  cusp  of  u,  has  touched  x  =  D  ve  obtain  figure  7 
in  which  the  simple  wave  appears.   We  have  only  drawn  the  lines 
with  slope  of  largest  magnitude  on  which  the  simple  wave  is  essen- 
tially constant.   All  tangents  to  the  slow  locus  of  U2  with  slope 
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smaller  in  absolute  value  than  that  of  these  lines  also  carry 
simple  wave  disturbances.   Pigires  8  and  9  show  the  disappearance 
of  the  cusped  region  of  the  slow  converted  to  fast  locus  and  the 
speeding  up  of  this  locus  and  the  slowing  down  of  the  fast  con- 
verted to  slow  locus.   In  figure  10  the  fast  converted  to  slow 
locus  and  slow  converted  to  fast  locus  have  crossed  and  the  former 
steepens  at  y  =  0 .   Finally  a  cusped  region  appears  in  the  fast 
converted  to  slow  locus,  figure  11.   In  the  limit  of  very  large 
time,  figure  12,  the  loci  approach  limiting  positions  with  some 
curves  approaching  each  other. 

With  some  knowledge  of  the  geometry  of  the  wave  fronts,  we 
examine  some  analytic  properties  of  the  solution  we  have  found. 
We  start  by  examining  the  wave  fronts  a  little  more  carefully. 
We  recall  that  for  t  fixed  (31)  gives  a  parametric  representation 
of  a  wave  front  with  parameter  p  where  h(x,y,t,p)  is  given  by  one 
of  the  forms  (2?)  through  (30).   We  readily  find  that  on  the  wave 
fronts 

ax  a^  2^    _      J       ^S^^'P^    ,   ,    ^\(i>P^  ]    /^^,•(l>P) 

3^=    -^/II3F-    -|^  3p2  ^^o  3p2  |/  ap 


>(32) 


ap  -      ap  ax/  ay        ^^j^I'P^  y  ^^2      "  ^o        ,^2        j/        ap        ^ 

where  the  indices  j  and  k  are  equal  to  3  or  4  depending  on  the 
forms  (27)  —  (30).   The  representation  is  singular,  corresponding 

.2        g     2 

to  possible  singular  points  of  the  curves,  when  (g— )   +  (-^j  =  0 ,  or 

.2       a^s  (l,p)       a^s  (l,p) 

^  =  X  ^ +  x  ^-5 =  0    .      (33) 

~  2  2         O       2 

ap  3P  8P 

Equations  (33)  and  (31)  determine  the  position  of  all  the  cusps 
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found  on  the  wave  fronts . 

Now,  let  x,y  be  a  point  on  a  wave  front  corresponding  to  the 
parameter  value  p.   Let  us  take  a  neighboring  point  x,y  where  x-x 
is  small  and  see  if  h(x,y,t,p)  has  any  zeroes  near  p.   On  ex- 
panding we  find 

0  =  s.(p)(x-x)  .  ^  IPzill  .  |!h_  (,_5)(p_p)  .... 

so  that  to  leading  order  and  by  use  of  (32)  we  find 


/+2s(l,p)  ,      -. 

(P-P)  =M  3s^l,p)  ^^Q^^-^)    • 

-^ ^P 

9s  (l,p) 

Now  for  j  =  3  or  4  and  p  positive  imaginary,  s.(l,p)/ — ^^^^ is 

J  op 

positive  imaginary.   Hence,  if  x(p)  is  increasing  as  Im  (p)  in- 
creases, then  the  root  of  h(x,y,t,p)  near  p  is  pure  imaginary  for 
X  >  X  and  has  non— zero  real  part  for  x  <  x.   Thus,  we  infer  that 
roots  of  h(x,y,t,p)  switch  from  real  to  imaginary  or  conversely 
across  a  wave  front  depending  on  whether  x  is  increasing  or  de- 
creasing as  the  parameter  p  increases.   Finally,  for  |y|  suffi- 
ciently large  we  infer  that  all  roots  of  h(x,y,t,p)  =  0  require 
p  pure  imaginary.   With  this  information  we  can  readily  track  the 
roots  of  h(x,y,t,p)  =  0  for  the  various  waves. 

To  illustrate  the  process  of  tracking  roots,  let  us  take  the 
simplest  case  of  u, ( | x-x  |,y,t).   A  root  can  only  switch  from 
pure  imaginary  to  complex  as  x,y,t  moves  across  a  wave  front. 
Since  for  |y|  sufficiently  large  the  roots  are  imaginary,  they  are 
Imaginary  outside  the  fast  locus  and  u,  i  0.   As  we  cross  the  fast 
locus,  one  complex  root  with  Re  p  >  0  appears  and  u,  ?  0.   On 
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crossing  the  slow  locus,  we  see  that  a  complex  root  disappears. 
Since  there  Is  only  one  complex  root,  we  conclude  that  inside  the 
slow  locus  there  are  only  imaginary  roots  and  again  u.^^  h  0.   Thus, 
we  have  shown  the  presence  of  a  lacuna  inside  the  slow  locus. 
Exactly  the  same  reasoning  shows  the  presence  of  a  lacuna  inside 
that  part  of  the  slow  locus  of  Up  not  covered  by  the  simple  wave. 

We  see  a  more  complex  picture  when  we  consider  u-.  corresponding 
to  the  fast  converted  to  slow  and  slow  converted  to  fast  loci.   In 
figures  6,  7,  and  8  there  is  a  lacuna  in  u^  inside  that  part  of  the 
slow  converted  to  fast  locus  that  has  two  or  three  cusps.   As  these 
cusps  disappear,  so  does  the  lacuna.   Some  of  the  time  the  lacuna 
in  U-,  overlaps  the  lacunae  of  Up  and  u-j^  and  corresponds  to  a 
lacuna  of  the  whole  solution.   At  other  times,  the  lacuna  is  in 
one  part  of  the  solution  only.   No  lacunae  are  then  present  in  u^ 
until  the  fast  converted  to  slow  locus  develops  a  cusp,  figure  11 
and  figure  12,  when  again  a  lacuna  reappears.   For  sufficiently 
large  time  there  is  a  lacuna  for  the  whole  solution;  at  other 
times  parts  of  the  solution  have  partially  non— overlapping  lacunae. 

The  analysis  of  the  behavior  of  the  function  h(x,y,t,p)  also 
indicates  that  as  one  approaches  a  wave  front  from  the  side  on 
which  it  is  convex,  the  fundamental  solution  tends  to  infinity. 
Thus,  much  of  the  disturbance  is  concentrated  near  the  wave  fronts. 

The  last  point  we  treat  is  the  behavior  of  the  solution  near 
the  very  singular  cusped  point.   We  wish  to  determine  how  singular 
the  fundamental  solution  is  near  the  cusp  and  whether  a  delta 
function  is  present  for  a  problem  in  magnetohydrodynamics .   Since 
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this  is  an  essentially  local  analysis,  we  shall  expand  functions 
as  necessary  and  take  only  the  leading  order  and  the  most  singular 
contributions.   We  will  not  carry  the  analysis  in  full  rigor,  but 
it  should  be  clear  how  to  fill  in  the  gaps.   We  take  x,y,t  in  the 
neighborhood  of  the  cusp  and  we  expand  for  p  small,  the  values  of  p 
that  determine  the  behavior  near  the  cusp, and  we  introduce  the 
function 

a+i<» 


k(x,y,t)  =  Re 


2tt 


^^K^'   -      A 


°  i-o^^o  ^PS 


dp 


O-loo 


2  .  2  c  „  /  o   o 
a  A  l-2p/  „  ^ 
o  o  ^  ^■•1     2.2 

a  A 
o  o 


a'^+A'^ 
o   o , 


)(y.   Et   _£X   /^_0(^, 


2.2 
a  A 
o  o 


2(a^+A^) 
o   o 


)) 


(3M 


so  that  by  reference  to  (19)  and  (26)  we  infer  that  near  the  cusps 


u-]_(x-x^,y,t)  =  k(x-x^,y,t)  + 


2.  „2 
2^2   2  5(y)H(t-(x-x^)  j-f^]    + 


a^+A^  „ 
o   of 


/a^+A^ 
/  o   c 


a  A 
o  o 


a  A^ 
o  o 


(35a) 


and 


U2(x+x^,y,t)  =  k(x+x^,y,t)  +  J -j-^   |  6  (y  )H  (t- (x+x^  )  j-§-^)  +  •• 

a  A  3-^ ^« 

o  o  o  o 


(35b) 


where  ...  implies  less  singular  terms.   Thus,  we  concentrate  on 
k(x,y,t).   We  introduce  new  variables  w,  X,  m  by 


p  =  ICJ 


X  =  t 


2a2 
a  A 
o  o 

2   2 

o   o 


,2fl2 


/a^A^     ,     a^A^   2 
-o-^^o         ^o^^o 


a^A^   2 
c  o  O  -1 
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where  X   and  m  are  the  rescaled  x  and  y  and  we  take  X   and  u  small 
to  be  near  the  cusp  and  we  find 

(36) 


k(x,y,t)  =  Re ^ 


dw  — 

-°°-10  U) 


o  ■  O 

We  recall  that  we  wish  to  pick  up  contributions  of  poles  in  the 

lower  half— plane  only.   We  see  in  (36)  the  origin  of  the  complex- 
ity:  M,  A  and  w  are  simultaneously  small.   We  can  regularize  the 
representation  of  the  roots  of  the  cubic  if  we  introduce  the  new 
variables 

Y 


:-M  +  Ju2-x3) 


1/3 


1 


6  =  (-P  -  y^) 


1/3 


(37) 


where  y   and  6  are  both  real.   Since  k(x,y,t)  =  0  if  the  cubic  has 

only  real  roots,  we  restrict  our  attention  to  the  case  of  complex 

2     ^ 
roots,  y  _>  X    ,    so  that  y   and  6  are  both  well— defined  and  real. 


The  roots  of 


■^ —  -  ■'h: —  +  U  =  0 


are 


Y+6,      -  \iy+6)    +  i|2(Y-6), 
Since  y-<5  >  0,  the  root  we  wish  is 


r(Y  +  6)  - 


./3 


(y-<5) 


V3/ 


=  -  ^(y  +  6)  -  ^{y-S) 


(38) 
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Thus,    we    obtain 


k(x,y,t)    =   Re 


1       (-0      o^2  o      o  0*= 


(A-^l 


Sttx    ^a   A 
o   o 


A      /a+A^   (w-^-Xu)) 


O     O 


where  w  is  given  by  (38).   On  expressing  k(x,y,t)  in  terms  of 
Y  and  6  through  (37)  and  (38)  we  find 


k(x,y,t)  = 


-1 
/Sttx 


(a^+A^)  2 
[   °   °  ] 

o  o 


1     MaX)^^f-^g^o(^^^^^)/2     2  ,3 


2  .2, 


2.2. 


a  A/A^a2  ^^y   -'    )-^3Y  6  ](y-6) 
o  ov  o   o 


,  u2<x3 


(39) 

2     3 

If  y   f.  ^   then  we  are  inside  one  of  the  cusped  regions  where 

2   3 
k  =  0.   Equivalent  to  the  restriction  y  -X^    >_  0  is  the  condition 

cited  above 

Y  -  6  >_  0 

The  singularity  in  k(x,y,t)  when  y  =  iS  is  the  singularity  on  the 
wave  front.   The  singularity  at  the  tip  of  the  cusp,  which  is 
additional  to  the  singularity  associated  with  being  on  a  wave 
front,  manifests  itself  In  k(x,y,t)  being  a  sum  of  terms  homoge- 
neous of  degrees  -3  and  -2. 

In  order  to  Interpret  (39)  we  take  a  convolution  of  k(x,y,t) 
with  data  appropriate  to  a  magnetohydrodynamics  problem.   That  is, 
we  Interpret  f  as  the  x  derivative  of  a  function  and  g  as  a  linear 
combination  of  x  and  g  derivatives  of  a  function.   Rather  than 
use  the  variables  x  and  y,  it  is  convenient  to  use  the  variables 
Y  and  e  =  &/y.      To  leading  order  In  small  terms 
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2   2.2 
o  o 


2d.  _  A     rt 
r  0  o  1 2 


d _  d      c  O  O  1  . 


O    O 


and 


3 


d   _  2  c  o  o  1 2 

9y    3x*-  2_^.2''  T7 

■^      a  +A  Y  ( 

00  ' 


2  9   ^  „ 


r;^!"'  '^ '  ^'^  ^^ 


so  that  If  we  combine  all  the  convolutions  into  one  form  of 
interest,  we  introduce 


J  = 


dx'dy'  k(x',y',t)  I    A  -^  +  B  — 

L   3x      9y 


i>i^   ,y   )  , 


or 


J  =  - 


/3 


^2 
o 


dY 


de  - 


dY 


de 


0    -00 


(a^+A^)f£+ga^Y(l+e^)/2 
Y^(l-e)(l  +  e^  +  e'^) 


2TTa  A  ,/a^+A' 
o  ov  o   c 

.  (A(Y^(l  +  e3)|j  -  Y^e^  |^)+B(2y£  ^  -   Y^  |y) ) X  (v^e ,y^ (l  +  e^) )    (40) 


where 
x[Y^t,Y^(l+e3)] 


\ 


a  A       "^  a  +A 
00         00 


-t 


/a^+A^   a  A 
o   o  r  o 


o  a  A    2  „ 
3  r  o  o  1  .2. 


1     V  2.2 
-1     ^  a  A 
o  o 


k  2   2 
a  +A 
o   o 


o_>2  Y^(l+£^) 
?J       3 


T  a  A     T 
1  -  |(-|-^)  yh 


so  that  Y  =  0  is  the  position  of  the  cusp.   We  assume  that  our 
data  is  non— zero  only  near  the  cusp  so  that  for  |y|  or  |e|  large 
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the  functions  vanish  identically. 

If  we  wish  to  recover  from  (40)  the  solution  to  the  initial 
value  problem  where  data  on  f(x,y)  is  given  we  set  f  =  1,  g  =  0, 
A  =  1,  B  =  0  and  :r^  =  f(x,y)  .   Similarly,  if  we  give  data  on 
g(x,y)  =  ^  +  ^  we  may  construct  it  as  the  sum  of  two  expressions 
of  the  form  of  J,  the  first  with  f=0,  g=l,  A=l,  B=0  and 

i(j  replaced  by  ()),  and  the  second  with  f  =  0,  g  =  l,  A  =  0,  B=l 

•x, 
and  \p   replaced  by  <}> . 

We  note  that  the  singularities  in  the  integrand  as  y  "^  0  or 

as  e  -»■  1  are  only  apparent,  as  can  be  seen  by  carrying  out  the 

1   .   -d 


implicit  differentiations.   If  we  replace  y—  by  , 


log  ll-e I , 

then  we  may  integrate  by  parts  with  respect  to  e  and  no  end  point 
contributions  appear.   Certain  of  the  remaining  terms  are  then 
exact  derivatives  with  respect  to  y  and  we  may  integrate  them  to 
find 


J  =  - 


/3^(t/^-^,0) 

^a  A 
o  o 


2iTa-A 
o  o 


H 


+A 


\    (a?+A5)f 


o   o 


de  log  I l-£ 


d_ 
de 


■    2  ^i 
1+e  +e 


1 


]    +   -^       de  log  |l-e|  ^  [ 2 V     ( 

'^      J  ^      1+e  +£    J 

-_oo  / 


dyde 


I    3x      3y       ^„2       9x9y       ^ 


'± 


(41) 


\  "-  "0^  3x"  """>^  8y' 
where  the  convolution  kernels  are  all  absolutely  integrable  and 
all  go  to  zero  as  y  tends  to  zero.  In  (4l)  we  have  isolated  in 
the  first  bracket  all  the  delta  function  behavior.   We  may  readily 
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evaluate  the  integrals.   Thus, 


de    log  |1-£|  ^  ( 


1  +  e  +e 

—  CO 

I 

and  if  we  set  e  =  1/e  we  find 

.3 


de  s."2rTrJ  -  P-^'    t::?  >- 


2 5"^ 

1+E  +e 


P.V. 
and  next, 
de  log  I 1-e 


de 


QE    C 

1-e  *- 


1+e 


d_ 
de 


1+e  +e 


e3 
2T~f 


de 


1  +  e- 


^)  =  -  p  V    ^^  r  ^^^'  ,  )  =  0 


1+e  +e 


P.V. 


1+e  +e 


de 


ae      (        e    ^ 

1^  >■ 2 — ?^ 

-^  ^   1+e  +e 


=  Re 


de 

1-e 


Im  e>0 


e3 


/3 


Hence,  the  delta  function  component  of  J,  Jr,  is  given  by 

4/[t 


^6  =  - 


f  I'^X 

2  ^aV 
o  o 


a^A^ 
^o^o 


and  by  reference  to  (35a)  and  (35b),  we  see  that  no  delta  function 
is  present  in  the  cusps  for  any  part  of  the  fundamental  solution. 

Finally,  it  is  of  some  interest  to  obtain  an  estimate  for  the 
kernel  k(x,y,t)  in  (39)  in  the  simpler  X,u  variables,  which  are 
essentially  Cartesian  coordinates  with  origin  at  the  cusp.   We 
Introduce 


k  = 


Cy5  -  D^^^-^^^V2 
[(y^-62)^+3y^62](y-6) 


(42) 
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^^^      1 ^       y^+y6  +  6^ 

and  we  may  determine  positive  constants  C-,  and  Cp  so  that 

'^  Y  +Y  5  +6 

Hence,  we  may  write 

CX  f,(X,y)  +  Dy  f2(X,y) 
k  =  , 


y,2.,3(2,2_,3)l/3 


where  f-,(A,y)  and  fp(X,)j)  are  both  positive  functions  uniformly 
bounded  away  from  zero  and  infinity.   We  conclude  that  the  nature 

of  the  singularity  near  the  cusp  is  given  essentially  by  the 

^   ^-      ^        fo   ^    ,3n-1/3   ^    p   ,„  2  ,3n-1/3 
functions  (2y  -X    )  and — (2u  -^    ) 


17^  JT^x 


3 


Although  the  two  functions  look  similar,  we  recall  that  the  deriv- 
ative of  the  first  with  respect  to  A  contains  a  delta  function 
while  the  A  and  y  derivatives  of  the  second  have  none.   We  see 
that  the  behavior  near  A  =  y  =  0  is  strongly  non— uniform.   One  can 
easily  plot  the  level  lines  of  the  functions,  but  the  pictures  are 
not  very  informative. 
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fig. I  Branch  line  and  connections  of  the  functions 
Sj(k,p)  in  the  complex  p  plane. 


Branch  Point 
Branch  Line 

Curve  on  which 
—  Res|(k,p)=0, 

i«i.2.3.4 


^ 


> 


Contour  CI 


Contour  C2 
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fig.  3     Wave  fronts  of  part  of  the  Fundamental  Solution 
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fig. 4     Wave   Fronts  at   t  =  l.00 
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fig.  6  Wave  Fronts  at  t=l.25 
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fig.  7    Wave  Fronts  at  t =1.3 3 
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